Coherent quantum oscillations in coupled traps with ultracold atoms 



N.M. Bogoliubov*, J. TimoncnT, and M.B. Zvonarev"^ 
* Steklov Istitute of Mathematics at St. Petersburg, Russia 
^ Department of Physics, University of Jyvdskyld, Finland 

The dynamics of two traps with ultracold atoms and connected by Josephson type coupling, is shown 
to exhibit a transition from dispersive dynamics to localized coherent oscillations. This transition 
is controlled by coupling strength and energy offset between the traps. The dynamics is also shown 
to be exactly that of a Heisenberg chain with a linear magnetic field. Possible applications for 
"quantum-state engineering" are discussed. 
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Trapping and cooling of atomic gases has produced an 
interesting laboratory for analyzing the transition from 
classical to quantum behavior. The fascinating new phe- 
nomena include the creation of Bose-Einstein conden- 
sates in magneto-optical traps, and the realization of cou- 
pled quantum systems with BEC's in [^,D, e.g., optical 
lattices. One can now analyze in some detail the proper- 
ties of Josephson currents [^Q and Bloch oscillations ||] 
that appear in the coupled systems, and entanglement 
of quantum states in individual traps. Possible use |^] 
of trapped ultracold atoms in realizing a quantum com- 
puter has also been considered (for "quantum-state 
engineering" with conventional Josephson junctions see, 
e.g., the recent review [^). 

It is evident that quantum dynamics of these systems 
has become ^ an important topic. This is particularly 
so if one tries to manipulate the dynamics as is needed, 
e.g., in the case of quantum computer. Also, rapid de- 
coherence of oscillatory dynamics even without external 
perturbations would make controlled operations very dif- 
ficult to perform. In coupled atom- field systems, such as, 
e.g., the one of two-level atoms interacting with a high-Q 
cavity, escaping photons form the main loss mecha- 
nism for quantum coherence independent of the internal 
dynamics of the system. It may thus be useful to consider 
systems in which controllable quantum dynamics can be 
realized without active involvement of a quantized elec- 
tromagnetic field. 

We consider in this Letter a coupled system of two 
traps with ultracold atoms, possibly in a condensed state 
so that they can each be described with one boson field 
6i_2,^i2; with subindeces denoting the trap. Experi- 
ments are now approaching the limit in which a fully 
quantum mechanical treatment of the problem is neces- 
sary, and it is actually this limit which is interesting from 
the standpoint explained above. First attempts have al- 
ready been made (l^Jl^ to analyze this system quan- 
tum mechanically, and we shall concentrate here on the 
quantum dynamics of the system, which so far has re- 
mained mostly unexplored. The two traps are assumed 
to be connected by a Josephson type of coupling, and 
atoms are assumed to interact via delta function poten- 
tials so that their interactions are j^jl^,|l^ of the form 
cb^jb^jbjbj] j = 1,2. Here c is the interaction strength 



which for analytical simplicity is taken to be of oppo- 
site sign for the two traps (this is now experimentally 
possible, and results are expected to be similar for gen- 
eral c albeit much more difficult to achieve analytically) . 
We have previously shown that for large Josephson 
coupling a natural quantum Hamiltonian for systems of 
this kind can be expressed in terms of the number op- 
erators Nj = b^jbj and the exponential phase operators 
f/,^. = {Nj + l)-i/2^^.^(^t ^ ^t(^^. ^ i)-i/2^ Notice that 

these phase operators are one-sided unitary: (l)j(j>j = 1, 

(j>'j(j>j = 1 — |0)j(0|j, with \0)j appropriate vacua. In this 
case of two coupled traps the Hamiltonian thus becomes 

Htr = cNl - cNl (1) 

+ {u; + A)Ni +ujN2-^ {^MI + 0l02) , 

in which lu and lo + A are the single particle energies 
(frequencies) in the two traps. We include an offset A 
in the energies ( "detuning" ) because this is an important 
parameter in the problem as will become evident below, g 
is the strength of Josephson coupling. The total number 
of particles is conserved, [Htr,Ni + N2] = 0, so we can 
consider without loss of generality the Hamiltonian Htr — 

7? = ATVi - I + 4^2) , (2) 

where A = A + 2c{Ni + N2). The ordinary (orthonormal) 
Fock states of the system can be formed from the vacuum 
states by operating with the phase operators, 

\N,,N2)^(4r'{4)'''mi\o)2 (3) 

such that {Ni + N2)\Ni,N2) = (A^i + N2)\Ni,N2) = 
N\Ni, N2). In problems like entanglement of quantum 
states or quantum computing one is interested in the dy- 
namics of a state which is not an eigenstate of the Hamil- 
tonian. The Fock states are not eigenstates of Hamilto- 
nian Eq. (H) , and in view of possible applications of the 
system, we consider next the dynamics of Fock states. 
We find it instructive to consider the dynamics in terms 
of the correlator 

^^i{t)^{j,N~j\e-^'"\l,N^l), (4) 
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which describes the time evolution of the occupation of 
particles in the two traps, and has an obvious connection 
to probability of observing a given state. 

This correlator can be solved exactly for the Hamilto- 
nian (j^). The most convenient way to find the solution 
is to first construct an "equation of motion" for the cor- 
relator, which for fixed subindex j takes the form 
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(5) 



and a similar equation is found for fixed I. The initial con- 
dition for the correlator is $^(0) = Sji, and the boundary 

conditions are = for j, I = -1, iV+1. Notice that 

this differential equation is a generalization of the equa- 
tion of motion for a lattice particle in a linear (electric) 
field 

Consider first the case of vanishing "field" or "effective 
detuning", i.e. A = 0. In this case Eq. (||) can easily be 
solved and we find 
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in which the (Bloch) spectrum Ek is given by Ek = 
— 5 cos(7rfc/(iV -t- 2)). The time evolution of any initial 
state can now be followed. We show in Fig. 1 the proba- 
bility for the initial state with equal number of 
particles in both traps, j = I = N/2. 
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It is evident from this figure that there is no coherent 
oscillation in the occupation numbers, the probability is 



rapidly dispersed into the available Fock space. The be- 
havior is analogous to the one found in |17| for a "spin 
wave" in an optical lattice, but there the dynamics was 
not followed very far (see also discussion below for con- 
nection to the Heisenberg chain). Notice that the time 
evolution is not translationally invariant at intermedi- 
ate time scales, and thus exhibits "aging" . There is 
however an almost complete recovery of the initial prob- 
ability distribution at a late time (which is repeated pe- 
riodically) . The small wiggles in the probability that are 
clearly visible at very early times are the "Rabi oscilla- 
tions" of |l5). 

For non-zero effective detuning A 7^ 0, Eq. (^ can 
also be solved exactly. For the same initial condition as 
above, 



fe=i 



(7) 



in which A = 5/A, and B{i', A; k) arc Lommel polynomi- 
als. In terms of the Bcssel functions of the first and sec- 
ond kind B(iy, A; fc) = ^[J,(A)r,+™(A)-J,+„(A)y,(A)]. 
The boundary conditions given below Eq. (0) mean that 
3(1/, A; A^-|-2) = 0, and the N+1 roots of this equation are 
the Vk- The spectrum Ek is now Ek = — A(j/fe + 1). For 
small N/X this spectrum is similar to the Bloch spectrum 
found above for zero effective detuning. There is a cross- 
over around N X into a linear spectrum: Ek ~ Ak for 
N/X ^ 1. This latter kind of spectrum is also called the 
Wannier-Stark ladder. 

For N/X <C 1 we thus find the probability |$^(t)P 
behaves very similarly to that for zero effective detuning 
shown in Fig. 1. Around X ~ N there is a transition to 
a completely different behavior. An excitation composed 
of a superposition of states that are localized in the Fock 
space is formed in the system such that coherent oscil- 
lation appears in the occupation probabilities of the two 
traps. 
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FIG. 2. Probability |<E>|?(t)|^ for j = 10 as a surface dia- 
gram in the tl plane for A = 5, and $ji'(0) = 5io,io- Height is 
indicated with a grey scale such that black means zero height. 



We show |$j;(t)r for N/X = 4 in Fig. 2, where these 
features of the solution are clearly visible. For increasing 
N/ A the number of localized states in the superposition 
decreases, and the probability becomes more and more 
localized around the center line (i.e. the initial value). 
This is evident from Fig. 3, which shows \^ji{t)\'^ for 
N/\ = 10. 
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FIG. 3. Probability |'l>|?(t)P 



for j = 10 as a surface dia- 



gram in the tl plane for A = 2. 



and $j?(0) = 5io,io- Height is 



indicated with a grey scale such that black means zero height. 

We could call the localized oscillatory solution as a 
"trapped quantum soliton". We believe that this solu- 
tion should be experimentally observable in coupled traps 
with ultracold atoms, in which the parameter A can fairly 
easily be controlled by controlling either A or g or both. 
Notice also that the states that form the superposition 
are entangled. This system would thus be suitable for an- 
alyzing entanglement, and perhaps even quantum com- 
puting as "engineering" of states is possible. 

It is also instructive to consider the limit in which 
one of the coupled traps becomes macroscopic (classical): 
iV — > oo (i.e. either A'^i or N2). For simplicity we consider 
this case only in the limit of weakly interacting particles 
and set c = 0. The correlator $°^(i) = (j| exp(-zti/oo)IO 
satisfies the same equation Eq. (Q) as for finite but now 
the boundary conditions are $°^(i) = for j,l — —1. 
The Hamiltonian becomes in this limit 



AiV~|(0-t-0t). 



(8) 



Notice that the (new) number operator N does not any 
more commute with the Hamiltonian (this is similar to 
the corresponding limit in the Bogoliubov theory p9| ] 



of Bose condensates). If we express this Hamiltonian 
Eq. (^) in the Fock space, we find that 



Hoo - A^n|n)(n| - f ^ {\n){n + 1| + |n+ l)(n|) , 



(9) 



This form of the Hamiltonian appears in many problems 
related to trapped gases of ultracold atoms, e.g., in the 
tight binding limit of a ring-shaped trap with a moving 
defect g. 

The correlator (t) can again be solved exactly with 
the result 
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where summation is over the Vk determined by the solu- 
tion of the "index equation" J^^. (A) = 0. The spectrum 
is now Ek + A = — A|i^fe|, which is very similar to the 
Wannier-Stark spectrum for small k. The solution is also 
in this case a coherent oscillation in the occupation num- 
ber of the finite trap similarly to the solution shown in 
Fig. 3. The same "quantum features" of the solution can 
thus be realized with only one trap with a small number 
of weakly interacting particles. 

For zero effective detuning coherent oscillations dis- 
appear as expected, and we find that asymptotically 
|$°^(t)| decays as t~^/'^. The dynamical exponent z = 
3/2 also appears in nonequilibrium growth processes with 
one-dimensional interfaces pC[ |. 

We shall finally make a connection to quantum com- 
puting even though a detailed analysis of this aspect of 
the problem cannot be given here. We show here that 
the correlators ^fi{t) considered above can be mapped 
exactly to the corresponding correlators of an isotropic 
Heisenberg spin 1/2 chain in a linear field, for which all 
logical operations of quantum computing have already 
been |2^ constructed. 

Consider the one spin-fiip correlator (related to the 
"NOT" operation) 
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(11) 



where the ground state |0) is that with all spins pointing 
up, and the Heisenberg Hamiltonian Hx:^^ is given by 



N 



<^n ^,t+l + ^(«+l - 1) 
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(12) 



Here tr^ , cr^ are Pauli spin matrices, and hn is an external 
mag netic field. The constant factors in the Hamiltonian 
( |l^ ) are added to ensure that Hxxx\^) = 0. As before we 
can derive an equation of motion for correlator (|l^) , and 
find for fixed index j that 
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-2{l + ^yf,+\t) (13) 

with a similar equation for fixed I. The initial condition 
to be apphed here is i^j)'^^(O) — Sji, and the boundary 
conditions appropriate for the lattice of + 1 sites are 
Fj^+^(t) = 0; j, Z = 0, + 2. As the constant factor in 
the last term on the right side of Eq. (^3|) only gives a 
constant global phase in the solution, it is evident that 
the absolute value of a solution of this equation for a con- 
stant magnetic field is identical with the corresponding 
solution of Eq. (H) with A = 0, and with that of Eq. ^ 
with A ^ for a linear magnetic field /i„ = nh. There is 
obviously a one-to-one mapping from the Heisenberg dy- 
namics into that of coupled traps with ultracold atoms. 
The number of lattice sites in the Heisenberg chain plays 
the role of total particle number in the traps. Thus, e.g., 
the limit A^ ^ cxd in the coupled traps, with Hamilto- 
nian (^, corresponds to the limit of semi-infinite chain 
in the Heisenberg problem. Notice that for a non-zero 
linear magnetic field, the localized coherent oscillations 
of Figs. 2 and 3 appear also in the Heisenberg chain. 

In conclusion, we have considered a model for two traps 
with ultracold atoms and with a strong Josephson cou- 
pling between them. We find that for non-zero effective 
detuning in the trap energies there is a transition to local- 
ized coherent oscillation in the occupation probabilities 
of the traps. Below the transition, i.e. for N/X ^ 1, 
the dynamics is dispersive as in the case of vanishing 
offset. We also found that there is a one-to-one map- 
ping between these dynamics and those of an isotropic 
Heisenberg chain in an appropriate magnetic field. This 
opens up the possibility in principle of realizing quantum 
computation, or simulation, in coupled-trap systems with 
only the detuning and the coupling of the traps as the 
relevant parameters. 
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